Strong quantum memory at resonant Fermi edges revealed by shot noise 
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Studies of non-equilibrium current fluctuations enable assessing correlations involved in quantum 
transport through nanoscale conductors. They provide additional information to the mean current 
on charge statistics and the presence of coherence, dissipation, disorder, or entanglement. Shot noise, 
being a temporal integral of the current autocorrelation function, reveals dynamical information. 
In particular, it detects presence of non-Markovian dynamics, i.e., memory, within open systems, 
which has been subject of many current theoretical studies. We report on low-temperature shot 
noise measurements of electronic transport through InAs quantum dots in the Fermi-edge singularity 
regime and show that it exhibits strong memory effects caused by quantum correlations between 
the dot and fermionic reservoirs. Our work, apart from addressing noise in archetypical strongly 
correlated system of prime interest, discloses generic quantum dynamical mechanism occurring at 
interacting resonant Fermi edges. 



Non-equilibrium electronic shot noise is a powerful 
diagnostic tool revealing properties of mesoscopic sys- 
tems inaccessible by the mean current measurements?" 3 -' 
For example, recent noise measurements in the Kondo 
impurities-^ have brought new insights into strongly cor- 
related transport. Shot noise is sensitive to the pres- 
ence of non-Markovian dynamics^ intensively studied in 
broad context ranging from photosynthesis to quantum 
informationPEl However, most theoretical proposals as 
well as the newest quantum-optical experimental studjff^l 
rely on extensive engineering and control of the system 
and/or environment (bath) and a clear observation and 
identification of the quantum memory effects in "natu- 
ral", i.e., routinely fabricated solid-state systems has not 
been reported yet. 

In resonant tunneling, which is ubiquitous in quan- 
tum electronic transport, the charge dynamics of the 
resonant level can be described by a simple Marko- 
vian master equatiorP as long as the relaxation time re- 
lated to the inverse of the transfer rates is long com- 
pared to the characteristic memory time of the fermionic 
bath (leads) given by the inverse temperature and/or 
the detuning of the level from the chemical potentials 
of the leads. Comparable time scales for system relax- 
ation and bath memory break down the conventional de- 
scription for low-temperature on-resonance transport and 
indicate^ strong non-Markovian features which, together 
with many body-interactions typical for small nanostruc- 
tures, influence the low-temperature width of the reso- 
nant steps in the current-voltage characteristics}-^ the 
decay of the level occupations,^ or the noiseP For the 
noise, significant deviations from the conventional mas- 
ter equation description have already been observed,^ 
although their origin has not been identified. 

In this work we present new experimental results on 
the low-temperature noise measurements in the Fermi- 



edge singularity (FES) regime^Hini together with a the- 
oretical analysis clearly revealing the presence of strong 
quantum memory around the edge. The noise-around- 
the-edge puzzle^ is briefly introduced in Fig. [l^, where 
the measured points are contrasted with the standard 
Markovian theory^ (black line) showing large deviations 
of ~ 15% in the Fano factor F = S/2eI (a convenient 
dimensionless measure of the shot noise S), far beyond 
the experimental uncertainty. Moreover, the measured 
dip breaks the Markovian lower bouncPof 1/2, which is 
a clear witness of strong memory. The blue line, nicely 
coinciding with the data, is our new theory accounting 
for the memory effects. 

The F ES, a paradigmatic exactly solvable many-body 
problem) 21 * 22 ' which originates from the Coulomb inter- 
action of conduction electrons with those on a localised 
discrete level represented by core shell electrons or quan- 
tum dot (QD) levels, was first predicted in the X-ray 
spectra of metalsP^ but its signatures are observed 
also in resonant tunneling set-ups as a (truncated) 
power-law singularity of the mean current / around, 
e.g., the emitter Fermi energyPZHSSI The interacting 
resonant level model describing the FES transport setup 
has served recently as an important benchmark for 
novel quantum transport techniques'^ including the 
noise calculation^ at its exactly solvable self-dual point 
(different from our experimental regime). 



Results 

We first describe cross-correlation measurements of cur- 
rent shot noise in self-assembled InAs QDs in the FES 
regime. The experimental set-up is depicted in Fig. [TJi 
and explained in more detail in the Methods section. At 
zero bias voltage the ground state energy level of the 
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FIG. 1: Shot noise power measurement around the Fermi edge singularity, (a) Simplified schematic of the studied 
device consisting of an InAs QD (pyramid) between emitter (E, blue) and collector (C, red) and the equivalent detection circuit, 
(b) Energy levels of the dot and leads. Left: zero applied bias with unoccupied dot level high above the lead chemical potential. 
Right: threshold bias when the dot level aligns with the emitter chemical potential and the current strongly enhanced by the 
Fermi edge singularity starts flowing, (c) Mean current / (solid line) and shot noise power S (symbols) as functions of applied 
voltage at T = 70 mK and B = 0T. (d) Energy dependence of ie and 7c in the off-resonant regime as determined from the 
measured current and shot noise (symbols outside the shaded region) and around the resonance with 7c fixed and 7_b calculated 
from the mean current (lines within the shaded region) . (e) Fano factor around the resonant edge (corresponding to the shaded 
range of (c) and (d)). Experimental values supplemented with their estimated errors (details in the Methods section) are 
contrasted with the Markovian approximation based on tunnel rate values from (d) (black line) and the full non-Markovian 
theory (blue line). 



InAs dots lies far above the emitter Fermi energy [Ie (see 
Fig. [I]d, left). Therefore, a large threshold voltage bias 
Vth « 170 mV applied to the collector lead is required to 
shift Ed to resonance with the emitter Fermi energy by 
electrostatic gating with the leverage factor r\ = ^yZy" \ 
giving the fraction of the bias voltage dropped at the 
emitter-dot junction (see Fig. [lj>, right). On resonance, 
the tunneling current sets in and displays a sharp peak 
(shaded part of Fig. solid line) due to the Fermi edge 
singularity caused by the Coulomb interaction of the oc- 
cupied dot level with the electrons in the emitter lead 
(there is no relevant interaction with the collector due to 
the asymmetry of the setup). Further increase of the bias 
causes a decrease of the current due to the decrease of 
the emitter rate induced by the three-dimensional den- 
sity of states (DOS) in the emitter!^ Together with the 
current-voltage characteristics on a large voltage scale, 
i.e., far around the edge, Fig. [It shows the measured 
shot noise power S (symbols). 

Far enough from the edge, i.e., outside of the shaded 
region of Fig. [lb, we can use the standard Markovian 



master equation and evaluate the emitter je and col- 
lector 7c tunneling rates, Fig. [TJi, from formulas^ I = 
2e^Elc/{^lE+lc) and F = 1- ^eIg I '( 2 7£+7c) 2 (ex- 
cluded double occupancy due to strong onsite Coulomb 
interaction implies usage of 2^/e instead of just je as 
for a noninteracting resonant leveP^). While the collec- 
tor rate jc = 8 • 10 9 s _1 is basically constant, te re- 
flects the energy dependence of the emitter DO£p2l and 
exhibits an expected asymmetry of the tunneling barri- 
ers with 7b/7c ranging between 0.06 and 0.22. Plausibly 
assuming constant 7c throughout the resonance we can 
analogously to Ref. [TH] extrapolate the -fE to the reso- 
nance (shaded) region in Fig. [TJi (solid lines) from the 
expression for the current. The resulting Markovian pre- 
diction based on these extrapolated rates (black curve in 
Fig. [l^) clearly exhibits substantial deviations from the 
measurement inexplicable by experimental errors. Obvi- 
ously, the resonant transport regime calls for a radically 
new theoretical understanding. 

Discussion 

Using the procedure briefly described in the Methods sec- 
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tion for B — 0, we arrive at a non-Markovian generalised 
master equation (GME) for the occupations of the reso- 
nant level, where pi(t) is the probability that the level is 
occupied by an electron, while po(t) = 1 — Pi(t) denotes 
the probability of the dot being empty, 



dpo(t) ,,s 
__ =7oRl ( t ) + 



dt'[ 7 Ut')pi(t-t')-H(t')po(t-t% 

The expressions (K = 1) for the forward/backward non- 
Markovian electron transfer rates across the QD/emitter- 
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lead interface 

volve standard FES Green's functions, whos e eval- 
uation is a known result of the FES theory^ 22 * 30 1 
This results in the explicit form of the rates enter- 
in the Laplace space they read 

' R ( l—QL 1 z+7c(l+'A)/2 
D \ 2 ' 2nk B T 

1e {z; -A), where A = = 

6T '^7c/2 th ' > * s t ^ ie dimensionless energy/voltage distance 
from the resonant edge, a is the FES critical exponent, 
and B(x,y) denotes the beta-function. In the zero- 
temperature limit the formula simplifies to j E (z; A) oc 

z + 7c (i+iA)/2 J I ■ When the counting field x at 
the emitter junction is includecP 3 -^ the GME memory 
kernel corresponding to equation ([I]) is of the form 



_ ( -2 1e {z; A) 7 |(z;A) e -x+ 7c 



2 7 £(z;A)e* - 1e (z;A)- 1c 



(2) 



Using the standard procedure for the cumulant 
evaluatio n^ 31 ! on this memory kernel, we get, using 
the abbreviations r y E b = 7^/ 6 (0;A) and j E b ' = 



d 1 f E /b (z;A)/dz\ 



= o * the formulas for the mean current 
and for the non-Markovian Fano factor 



F = 1-4- 
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The last term in the Fano factor, proportional to the 
derivatives, constitutes the non-Markovian correction. 
Well above the edge, where both the back-flow j E and 
the non-Markovian features can be neglected, we recover 
the standard master equation result (with 7^ = 7^). 



In the lowest order in r y E ' b /jc we can write F ~ 
1 - 4 7^/7c + 47^, with magnitude \r E h E \ « 
[max(fcsr, 7c, 7cA/2)] . This implies that for low 
temperatures ksT < 7c and close to the edge |A| < 2 
the non-Markovian correction is governed by the collector 
rate 7c and, thus, it is of the same order as the Marko- 
vian correction —^i E j^ic to the Poissonian noise (with 
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FIG. 2: Mean current fits for FES in magnetic field. 

Top (bottom) panel: I-V curves for different temperatures 
(colour-coded as shown by the middle panel and horizontally 
shifted with respect to each other for clarity) at B — 9.75 T 
(B — 0). Experimental data are plotted in colour and the 
black dashed curves are the theoretical fits for a unique set 
of 6 (4) parameters at all temperatures. Middle panel: Dif- 
ferential conductance dl/dV as a function of the bias voltage 
and magnetic field at the lowest temperature T = 70 mK. 
The Zeeman splitting of the edge is clearly visible. The top 
and bottom panels depict the cuts along the corresponding 
borders of the middle panel; the shaded stripes in the panels 
indicate matching ranges for the lowest temperature curves. 



F = 1) due to correlations caused by sequential occupy- 
ing and emptying of the QD. Being of quantum origin, 
it vanishes fast with increasing temperature kgT > 7c 
which kills quantum correlations between the dot and 
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FIG. 3: Memory effects on the FES noise. Top (bottom) panel: mean current (Fano factor) for two values of the 
magnetic field specified in the top panel and various temperatures shown in the bottom panel. Measured Fano factor with the 
estimated error-bars (explained in the Methods section) is compared to theoretical predictions based on parameters obtained 
from the fits of Fig. [2] Insets: Comparison of non-Markovian theory (solid lines) with the Markovian approximation (dashed 
lines) for corresponding magnetic fields (individual insets) and temperatures (curves within insets; horizontally displaced for 
clarity). Differences between the two curves are highlighted by colours according to their sign. Left detail: zoom onto the 
low-temperature Fano factor curve in the region around the upper Fermi edge shown by the dashed rectangle. Measured data 
with their error-bars are supplemented with both the non-Markovian as well as Markovian predictions in the spirit of the insets. 



leads responsible for the memory effects. Moreover, it 
generically assumes both signs — negative above the 
edge, further suppressing the Fano factor as in Fig. [Tjj, 
but also positive below the edge in the purely quan- 
tum tunneling regime, where it counteracts the classical 
term by increasing the Fano factor to potentially super- 
Poissonian values (F > 1). While the noise suppression 
can be achieved by memory of any origin, quantum or 
classical, the noise enhancement is a fingerprint of subtle 
quantum correlations. 

We now demonstrate these concepts by more elabo- 
rate analysis of the experimental data acquired at var- 
ious values of the temperature and magnetic field. We 
start by fitting the experimental data for the mean cur- 
rent (insensitive to memory) around the edge(s) with a 
straightforward extension of the above theory to the case 
of two spin-split levels due to the magnetic-field with re- 



sulting 3x3 (double occupancy excluded) memory kernel 
analogous to equation ^ as shown in Fig. [2] In the 
finite magnetic field case 6 free parameters were fixed 
by fitting simultaneously curves at various temperatures, 
namely two independent critical exponents = 0.40, 
a± = 0.43 and thresholds Vtht,^ together with an overall 
prefactor to the emitter rates and the leverage factor 77, 
while in the B = case only 4 parameters due to a single 
resonance peak were fixed with a = 0.28 in qualitative 
agreement with in-depth investigationsP2-_0 

By this procedure all parameter values are fully de- 
termined and the predicted Fano factor curves in the 
lower panel of Fig. [3] are free of any ambiguity. Consider- 
ing this, the correspondence between the measurements 
(points with error-bars) and our non-Markovian theory 
(lines) is quite remarkable in all cases encompassing two 
magnetic field values and various temperatures. We also 
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compare the Markovian, i.e., with the derivative terms in 
equation ([3]) omitted (dashed lines), and non-Markovian 
(solid lines) predictions in the insets and the detail of 
Fig. [3] with clear demonstration of the already mentioned 
non-Markovian features in the low-temperature Fano fac- 
tor, namely, the significantly more pronounced dip on the 
high-voltage side of the FES and the potentially super- 
Poissonian peak on the low-voltage side with fast destruc- 
tion of the non-Markovian corrections with temperature 
or distance from the resonant edge. 

All these features are clearly seen in the experimental 
data as well. The super-Poissonian Fano factor due to 
quantum coherence at the lower edge is not reliably 
confirmed experimentally because of associated large 
errors resulting from a ratio of very small values of both 
the current and noise (tunneling regime). Nevertheless, 
the experimentally observed peak just below the upper 
edge (see the detail in Fig. [3]) , although sub-Poissonian, 
is caused by the very same mechanism and is thus an 
indirect confirmation of the purely quantum memory 
effect. Altogether, the importance of the non-Markovian 
corrections due to quantum memory is established both 
qualitatively and quantitatively. 

Methods 

Experimental details. The studied InAs QDs are embed- 
ded in a GaAs-AlAs-GaAs resonant tunneling device pat- 
terned into pillars with a cross section of 9 x 9 /im suffi- 
ciently small to resolve single dot tunneling!^ The effective 
AlAs barrier widths of 4 and 3nm are slightly asymmetric. 
The measurements, whose schematic of the electronic setup 
is shown in Fig. [l^, were performed in a dilution refrigera- 
tor at temperatures down to 70 mK and magnetic fields up 
to 13 T. The DC-part of the source drain current I is mea- 
sured with a transimpedance amplifier which also biases the 
sample. Two 4.7 k£7 resistors convert the fluctuating current 
to voltages which are measured in a cross-correlation config- 
uration. Together with parasitic capacitances these resistors 
form RC-circuits which define the bandwidth of our experi- 
ment. To increase this bandwidth we use home-built coax- 
ial cables thereby lowering the total parasitic capacitance to 
20 pF. 

The voltage fluctuations are amplified by a two-stage low 
temperature amplifier based on the ATF34143 HEMT with 
a gain of 22 dB, followed by a room temperature amplifier 
with a gain of 60 dB. The amplified signal is filtered and digi- 
tised, and from the Fourier spectra the cross-correlation noise 
power Sab is calculated and averaged over 8 minutes. To 
retrieve the shot noise power the real part of Sab in the fre- 
quency range from 500 kHz to 3 MHz is evaluated. A technical 
background-noise, largely dominated by thermal noise sources 
like the first transistor stage of the cryogenic amplifiers and 
the conversion resistors, is measured at zero current and sub- 
tracted. At finite sample impedances the partially correlated 
thermal background is estimated and also subtracted. The 
correlation gain parameters are determined by noise ther- 
mometry. The error bars in Fig. []J and Fig. [3] consist of 
the statistical error, the error of the estimated background 
and the calibration error. 

Theoretical details Hamiltonian of simplified spin-less 
model of a resonant level tunnel-coupled to two leads (emitter 



E and collector C) and Coulomb-coupled just to the emitter 
reads H = H QD +H E + H C +H T + VxdU with H QD = eod^d, 
H P = E fe|9 ekpclpCkf,, Hf,T = T, k/) (. t k0d f Ch e + t* kfi J kf} d) and 
y x = T, kE , k ' E Vk E ,k' E cl E c k , E , where /3 = E, C; d and c k: g 
are QD and lead annihilation operators, ed and e k/3 are the 
energies of the QD level and of the electrons in the leads, re- 
spectively, while t kfj describe the tunneling between the QD 
and the leads. The last term describes the scattering of emit- 
ter lead electrons on an electron in the QD and is responsible 
for the FES phenomenon. Since the Fermi level of the col- 
lector lead is far below the resonant level and 7c does not 
depend on energy close to the edge (Fig. [ljl) we can use the 
method by Gurvitz and Pragei^to exactly integrate out the 
collector lead. This leads to the equation of motion for the 
density operator a(t;n) of the dot and the emitter resolved 
with respect to the number n of passed electrons through the 
emitter/QD interface. After introducing the counting field y 
as a conjugate variable to n, one can write the equation of 
motion for a(t;x) = ^] n £r ft n ) e " X ! partly expressed in the 
block formp^l a = (cr o, en, C01, cio) T , 



dcr(t;x) 
dt 



= -i(H E {x)a(t; x) ~ o{t; X )H E {-xj) 








7c 
-7c 

-7c/2 








-7c/2/ 



(4) 



where H E (x) = H QD + H E + VxdU + E^C^e^Vc^ + 



-x/V 



d) 



H E (—X) i s the appropriately modified 



Hamiltonian of the emitter and QD including the counting 
field. 

The emitter lead can then be handled perturbatively in the 
tunnel coupling tk E following closely the derivation for dissi- 
pative double quantum dot from Ref. [3T] by first separating 
equation Q into four equations for the elements of a. Tracing 
out the electron states of the emitter lead in the two equa- 
tions for the evolution of the diagonal elements ajj (j = 0, 1), 
allows us to find the evolution equations for the generalised 
QD occupations Pj(t; y) = Tr E <jjj(t; y) 



dp /i(t;x) 



dt 



±7cPi(t;x) 
T 2e ±x/2 ]T Im [t k „ Tr E (c kE a 01 (t; X ))} 



(5) 



The equation governing the evolution of aoi(t; y) which enters 
Eqs. |5| contains the diagonal elements ajj(t; X ) of the total 
density matrix. In order to close the equations for pj (t; y) we 
perform physically motivated QD-state-resolved perturbative 



decoupling of the density matrix into ajj it; y) = pj (t; y) <8> of 

7E 



with gf = exp ( - [H E + jV x - ^ E E fej5 c{ E c kE ]/k B T) /Zf 
being the grand-canonical density matrix of the emitter lead 
at temperature T and chemical potential fi E when the QD 
is empty (j — 0) or occupied (j = 1). Thus, after mul- 
tiplying the forward rate by 2 due to the interplay of spin 
and Coulomb blockadep^ we find equation Q (reducing to 
equation for y = 0) with FES Green's functions reading 

Go(t) = E^^WSi^^^+^^e-*^*^], 
Gi(t) = £*„.*/, tk B t* k , Tr E \e 



,i(H E +V x )t 



c kK e 
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